Direct measurement of the ^(1020) leptonic branching ratio 
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The process e + e~ — ■» has been studied by SND detector at VEPP-2M e + e~ collider in 

the </>(1020)-resonance energy region. The measured effective <f) meson leptonic branching ratio: 
B(<t> -> l + r) = y/B(cf> -f e+e~) • B(cf> -> p+p~) = (2.89 ± 0.10 ± 0.06) ■ 10~ 4 agrees well with the 
PDG value B{4> — > e + e~) = (2.91 ± 0.07) • 10 -4 confirming p— e universality. Without additional 
assumption of p—e universality the branching ratio B(<f> — * /i + fi~) = (2.87 ± 0.20 ± 0.14) • 10 -4 was 
obtained. 



Truly neutral vector mesons play an important role in hadron physics due to their direct coupling to photons. 
This phenomenon is the basis of the phenomenological Vector Meson Dominance model which successfully describes 
electromagnetic interactions of hadrons. The key parameters of this model are V— 7 coupling constants. They can be 
extracted from the vector meson leptonic widths under the assumption that leptonic decay proceeds via one-photon 
annihilation of the quark-antiquark pair constituting the meson. Leptonic widths also determine the total production 
cross sections of vector mesons in e + e~ annihilation and are important for calculation of the hadronic contribution 
to the photon vacuum polarization jjj . 

The V-7 coupling constant is just one number per vector meson. Could these numbers tell us something non-trivial 
about the underlying QCD dynamics? Shortly after the 1974 "charm revolution" , Yennie noticed that independently 
of the vector meson flavor content the following relation holds (||| 

T(V -> e+e - )/ < e q > 2 sa 12 keV, (1) 

where < e q > is the mean electric charge of the valence quarks inside the vector meson V in the units of an electron 
charge. For p, to and <\> mesons this gives the famous rule: T(p — > e + e~) : T(u> — > e + e~) : r(</> — > e + e~) =9:1:2, 
which can be considered as an SU(3) symmetry prediction. The surprising fact here is a relatively high (~ 10%) 
precision of the 9:1:2 rule despite SU(3)-flavor symmetry breaking. Inclusion of charm gives even more badly broken 
SU(4) symmetry, but Yennie's relation remains valid with the same precision, which means that SU(4) symmetry still 
persists for the leptonic widths ratios! Inspired by this strange fact, Gounaris predicted T(T — » e + e~) = 1.2kcV Q] 
and was closer to reality than any other author J3|. Current experimental situation with leptonic widths || is shown 
in Table | 



TABLE I. The leptonic widths of vector mesons. 







< e q > 2 
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6.77 ± 0.32 


1/2 


13.5 ±0.6 




0.60 ± 0.02 


1/18 


10.8 ±0.4 




1.30 ±0.03 


1/9 


11.7 ±0.3 


J/il> 


5.26 ±0.37 


4/9 


11.8 ±0.8 


T 


1.32 ±0.05 


1/9 


11.9 ±0.5 



In the nonrclativistic potential model |3| the leptonic decay width is given by the Van Royen-Weisskopf formula 
§: T(V ->e+e-) = 16Tta 2 <e q > 2 |#(r = 0)\ 2 /M 2 . 

Equation (0) implies then that quarkonium wave function at the origin \I/(r = 0) is proportional to the meson mass 
My. Note that for Coulomb potential |$(r = 0)| 2 ~ M v , while the linear potential gives |*(r = 0)| 2 — My. So 
the leptonic widths tell us that the actual potential appears to be something in between. But even if we postulate 
such a potential, the relation (0) still has no simple explanation. For light quark systems like p, to, and <f> relativistic 
corrections are essential. There are also strong interaction corrections governed by the scale dependent a s . It was 
argued || that these corrections modify the Van Royen - Weisskopf formula in the following way: 
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F(V -f e+e") m 16ttc^ < e, > 2 |*(r = l/m,)| 2 (l - 0.36 a,(M v ))/M£. (2) 

Intuitively, appearance of the constituent quark Compton wavelength l/m q in (||) looks natural, because in relativistic 
theory a particle cannot be localized within a region smaller than its Compton wavelength ||. Thus we can expect 
the quark-antiquark pair to annihilate when approaching each other's relativistic extents ^]. But this intuitive clarity 
of (En doesn't make an explanation of the remarkable regularity of simpler, because (Q) shows that leptonic widths 
are sensitive to the both nonperturbative and perturbative aspects of QCD. Thus it is not surprising that the leptonic 
widths become a traditional touchstone for various quark models [p||Tc|]. 

This paper is devoted to the measurement of the leptonic branching ratio of the </>(1020) meson. There are two 
leptonic decays: cp — > e + e~ and — > fi + fi~. The fi—e universality implies for these decays that B(<j) — > (i + fi~) = 
B{4> — > e + e~) x 0.9993. Presently only the <fi — > decay was measured directly ( jll]]- [0)- There are two PDG 

values for this decay branching ratio ||. One of them B{<p — > /i + /i _ ) = (2.5±0.4) • 10~ 4 is based on the experiments on 
photoproduction of meson JTuO . Another value of the branching ratio B(<f> — * fi + (i~) = (3.7±0.5)T0~ 4 is obtained 
from e + e~ experiments [ p^|jl4|J17f . In addition the CMD-2 experiment |l^,[l6| has some preliminary results on this 
decay. One can see that the difference between two PDG values for the decay <f> — > is about 2 standard deviations 

and the accuracy of these results is relatively low. Current branching ratio B(cj> — > e + e~) = (2.91 ± 0.07) • 10~ 4 || 
is based on measurements of the 0-meson total production cross section in e + e~ collisions. It was obtained by 
summation of all 0-meson decay modes: <fi ~ ¥ K + K~, KsKl, 3ir, etc.. Up to now the accuracy of B(<f> — > e + e~) 
was much higher than that of B(<p — > /i + /i~). But there is a serious factor limiting the precision of B(<f> — * e + e~) 
obtained in such an indirect way. It is the interference between 4> meson and other vector states, which description is 
model dependent. Direct measurement of the 4> — > e + e~ decay in the e + e~ — > — > e + e~ reaction is difficult due to 
its small probability and huge background from the e + e~ — > e + e~ Bhabha scattering. 

The decay </> — > reveals itself as a wave-like interference pattern in the energy dependence of the e + e~ — > 

cross section in the region close to the 0-meson peak. The amplitude of the interference wave is proportional to 
B((j> — > Z + Z~) = y/ B{<j) — > e+e _ ) • 5(0 — > The accuracy of the B(^> — > Z + / _ ) measurement in this case is 

limited only by uncertainty in the calculation of the pure QED part of the e + e~ — > cross section. The 0.2% 

accuracy claimed in leads to 0.8% systematic error in the interference amplitude. Large statistics collected by 
SND detector in the vicinity of the cj> resonance allowed us to make direct measurement of the leptonic branching 
ratio B(<p — > l + l~) with the accuracy comparable with that of previous indirect measurements of B(cj> — > e + e~). 

Our previous study of the e + e~ — » cross section was done using the 1996 data sample with the total 

integrated luminosity of 2.6 pb _1 p"^ ]. In 1998 two experimental runs were carried out in the center of mass energy 
range E = 984 — 1060 MeV in 16 energy points. The collider operated with superconducting wiggler |Q allowing 
to increase the average luminosity by a factor of two. Higher luminosity led to relative reduction of the cosmic ray 
background. The total integrated luminosity AL — 8.6 pb^ 1 collected in 1998 corresponds to 13.2 • 10 6 produced <j> 
mesons. 

The SND experimental setup is described in detail in ref. p9[ . The main part of the SND is a spherical electro- 
magnetic calorimeter. The angles of charged particles are measured by two cylindrical drift chambers (DC). An outer 
muon system, consisting of streamer tubes and plastic scintillation counters, covers the detector. The integrated lu- 
minosity was measured using e + e~ — » e + e~ events selected in the same acceptance angle as the events of the process 
under study e + e~ — > fi + pT . The systematic uncertainty of the luminosity measurement is 2%, but its contribution 
to the systematic error of the interference amplitude estimated using the process e + e~ —> 77 is only 0.8%. 

The primary selection criteria for events were similar to those of our previous work Jl7|| : 

• total energy deposition in the calorimeter is more than 270 MeV; 

• there are two collinear charged tracks in an event with acollinearity angles in azimuth and polar directions: 

Aip \< 10°, I AO |< 25° and with the polar angles within 45° < 9 < 135° ; 

• event is not tagged as e + e~ — > e + e~ by e/ir separation procedure j2(J. 

To suppress the background from the processes e + e~ — > ir + ir~, 7r + 7r~7r°, K$Kl, K + K~ the outer muon system 
was used: a requirement for both charged particles to produce hits in muon system renders contribution from this 
background negligible. For example, the contribution from the process e + e~ — > 7r + 7r~ is about 0.2% in the 0-meson 
peak. 

The cosmic ray background was suppressed by restriction of the time r measured by outer scintillation counters with 
respect to the beam collision moment ]ri| : r |< 10 ns. About 30% of events selected by the cuts described above are 
still cosmic ray background. To determine the contribution of cosmic background more accurately the selected events 
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FIG. 1. The Aifi distribution in e e — > /i fi events 
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FIG. 2. The A(p distribution for cosmic ray events. 



were divided into two classes: (1) | A(p |< 5°; (2) | Aip |> 5°. The resolution in Aip is about 1°. The Aip distribution 
for e + e~ —> /i + fj,~ events (Fig. [I]) was obtained from the experimental data after strong cuts on a difference between 
time measurements by the muon system for both tracks. Almost all events belong to the first class. The second 

class contains only 1.7% of events. The Aip distribution for pure cosmic ray events collected in a special run 

without beams is shown in Fig. ^. The uniformity of this distribution is an artifact of our DC track reconstruction 
algorithm in which the origin of a charged track in the X-Y plane is fixed to the beam collision point. From Fig. |^ 
the ratio between numbers of cosmic ray events in the two classes was found k cs = jVf s /N^ = 1.028 ± 0.033. 

The number of cosmic ray background events in the first class was calculated for each energy point Ei by the 
following formula: Ni s (Ei) = k cs ■ T(Ei) ■ dN2/dT. Here T(Ei) is a data acquisition time for an energy point Ei, 
dN2/dT is the cosmic event rate in class two averaged over both 1998 experimental runs. The net number of 
events for each energy point was obtained by subtraction of the cosmic ray background: N^(Ei) = Ni(Ei) — N^ s (Ei). 
The errors of the numbers N'-^Ei) include the errors of Ni(Ei) and N£ s (Ei). 

Energy dependence of the detection cross section was fitted according to the following formula: 



'A (£) 



a™(E)=a (E).R(E) 
a (E) = 2Tra 2 p(E)(l - f3 2 (E)/3)/E 2 , (3) 

where a is the fine structure constant; (3(E) = (1 — 4 • m 2 ^/ E 2 ) 1 / 2 ; m^, T^, A,p(E) = m| — E 2 — iET(E) are the 
</>-meson mass, width and inverse propagator respectively; ctq(E) is Born cross section of the process e + e~ — » fi + fi~; 

= Qf Jr e l 'V>' — interference parameter. The modulus of the interference parameter is related to the leptonic branching 
ratio: = B{<p — > l + l~) ■ 3/a. The factor R(E) takes into account the detection efficiency and radiative corrections: 

R{E)=e» ^ {E) - 2 . (4) 



(TO (E) 



1 - z 



Here is the result of Monte Carlo integration of the differential cross section of the process e + e~ — > |2l| ] 
for our geometric cuts with the energy dependent probability for muons to hit outer scintillation counters taken into 
account. The uncertainty in the energy dependence of this probability adds 1.7% to the systematic error of Q^. The 
parameter e M represents energy independent factor in the detection efficiency. It is determined mainly by the cut on 
total energy deposition in the calorimeter. The value of — 0.84 ±0.01 was obtained using Monte Carlo simulation of 
the process e + e _ — > fi + fi^(^) in SND detector p2fl . but in the fitting procedure e M was left free. In the calculation of 
radiative corrections the interference parameter was assumed purely real and equal to Z = B(<f> — > e + e~) -3/a = 0.120 
with the PDG table value for B(<j) -»■ e+er). 
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FIG. 3. The measured cross section of the process e + e —> fi + /i . 
TABLE II. The results of the fit with ip^ = for two experimental runs. Only statistical errors are shown. 



Parameter 


PHL9801 


PHL9802 


Combined 


X J /NDF 
Qij,, HT 2 
£ M , HT 2 

B(cj>^i + r), io~ 4 


19.4/15 
12.1 ± 0.6 
83.1 ±0.3 
2.99 ±0.15 


11.3/15 
11.0 ±0.6 
82.5 ±0.3 
2.74 ±0.14 


33.8/30 
11.9 ±0.4 
83.1(82.8) ±0.3 
2.89 ±0.10 



The fitting was performed for two experimental runs independently. Fits with a free ip^ yield the interference phase, 
which is in good agreement with the expected zero value: (1) ip^ = (1.0 ± 2.8)°, (2) ip^ = (0.1 ± 2.8)°. Therefore 
the interference phase was fixed to ip^ = 0. The fit results presented in Table [0] show statistical agreement between 
two experimental runs. Therefore combined fit was performed to obtain the final results which are listed in the third 
column of the Table [fl| The values of £ M obtained in the fit and from Monte Carlo simulation agree well. The energy 
dependence of the measured cross section and the fitting curve are shown in Fig. ||. The systematic error of 
includes 1.7% from the uncertainty in the energy dependence of the probability for muons to hit the outer system, 
0.8% from the luminosity measurements and 0.8% from the calculation of the radiative corrections. The resulting 
systematic error is 2%. 

In conclusion we obtain the following cf> meson parameters from the measured value: B(<fr — » = (2.89 ± 

0.10±0.06)-10~ 4 ; B{cj> -»• e + e~)-B(4> -> a* + aO = (8.36±0.59±0.37) • 10~ 8 . This result is in good agreement with our 
previous one B{<j> -» l+l~) = (3.14 ± 0.22 ± 0.14) • 10" 4 @. Using PDG value of B(<f> -> e+e - ) = (2.91 ± 0.07) • 10~ 4 
we obtain: B{<p -> ^ + aO = (2.87 ± 0.20 ± 0.14) • 10~ 4 . The good agreement of B(<f> -> m + aO and B(<j> -> e+e") 
confirms the fx—e universality. 
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